On the Charge Asymmetry of the Like-Sign 
. Lepton Pairs Induced by B — B Production 

Asymmetry 



O 



A. Belkov, T. Ilitcheva, S. Shulga* 
JINR, Dubna 



> 

^ \ Abstract 

In Monte Carlo simulation of pp and pn interactions, it is shown that the charge 
C*~) ■ asymmetry of like-sign lepton pairs can be observed as a manifestation of B — B 

production asymmetry. In this way, the B — B production asymmetry could be 
studied experimentally without full reconstruction of i?-mesons. 

It is known that the asymmetry between decays of B and B mesons provides a signal 
for CP violation, but can also be caused by the production asymmetry between B and B 
mesons in pp or pn collisions. The physics origin of asymmetry of heavy-meson production 
in hadronic interactions is related to a simple effect of the valence quarks in the colliding 
o5 ' hadrons and non trivial dynamics of the hadronization process [HI2j- 



At the parton level, the b and b quarks are generated symmetrically since their pro- 
duction is described within perturbative QCD by diagrams of hard scattering of partons 
which always arise through the g —> bb or gg — > bb couplings. As the result of the strong 
interaction in the confinement region of QCD, the colored quarks produced perturba- 
tively in the hard scattering processes are transformed into colorless hadrons. In the case 
of b quarks, this transformation (quark hadronization) may introduce an asymmetry be- 
tween a B meson and its antiparticle if there is an asymmetry in the quark and antiquark 
flavours that are available in the remnants of the initial hadron for 5-meson formation 
(see Fig.H}. 

Quark hadronization includes non-perturbative QCD processes and is treated and 
simulated by PYTHIA package |3J in the frame of the string fragmentation model. The 
Monte Carlo studies by using PYTHIA (see Ref. j2]) show that the overall B — B produc- 
tion asymmetry is expected to be below one percent for pp collisions at LHC and becomes 
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Figure 1: The mechanism of the production asymmetry in pp collisions due to asymmetry in 
flavours of quarks and antiquarks with which the produced b and b quarks can combine to form 
B mesons 

larger at smaller energies, giving a few percent or even above ten percent for pp and pn 
interactions at HERA-B. Thus, an experimental study of the B — B production asym- 
metry should therefore be most feasible at HERA-B with proton beam of 920 GeV. The 
measurement of the asymmetry in the restricted phase space regions, i.e. as a function 
of transverse momentum, rapidity, pseudorapidity, or Feynman variable xp of B mesons, 
provides much more information on the non-perturbative dynamics of quark hadroniza- 
tion. According to Monte Carlo studies of Ref. the B — B production asymmetry is 
positive and of order 10% in the central region of phase space, but becomes negative and 
very large in the extreme forward/backward regions close to the remnants. 

In spite of the great importance of a better understanding of B — B asymmetry for 
CP-violation studies and as a test of hadronization models as well as a mean to constrain 
their parameters, no measurements of this effect in hadronic 5-meson production has yet 
been performed. The direct measurement of the B — B asymmetry requires a triggering 
and reconstruction of 5-decays. One of the possible ways to perform the asymmetry 
measurement is to use lepton pairs from doubly semileptonic decays B — > l + X and 
B — > l~X to trigger the S-meson signal. The main difficulties in this case are related 
with the high rate of background from leptons coming from decays of charmed particles, 
resonances, kaons and pions. A separation of e/i pairs avoids resonance decays, Drell Yan 
production and other processes which feed di-electron and di-muon channels. But even 
after applying the additional kinematical cuts the problem of background reduction is not 
completely solved because of the dominance of the e/x signal from the doubly semileptonic 
decays of charmed particles, mainly from decays — > l + X and — > l~X. 

In this paper we propose to trigger the 5-meson signal by selection of the like-sign 



lepton pairs 11. This makes it possible to largely suppress the background originating 
from decays of charm particles and, after applying an additional cut on lepton transverse 
momentum, the Z ± Z ± signal from doubly semileptonic decays of B mesons becomes clearly 
seen. Moreover, the B — B production asymmetry induces the charge asymmetry in like- 
sign lepton pairs which, therefore, reflects the non-perturbative dynamics of hadronization 
processes but can be measured directly without reconstruction of 5-meson decays. The 
detailed Monte Carlo studies of lepton pair production in pp interactions at the proton 
energy of HERA-B have been performed for this paper by using the version PYTHIA 
6.158. We analyze the contribution of various physical sources to like-sign lepton pair 
production and their asymmetries. A special study has been devoted to the role of B°—B 
oscillations which obscure the lepton-pair asymmetry induced by asymmetry of B — B 
production (see also discussions in Ref. [2j). 

In this paper we consider only pp interactions at fixed target with proton beam energy 
of 920 GeV. Generation of bb and cc events in PYTHIA is based on the description 
of heavy flavor production within the usual parton model. At HERA-B energy, this 
approach assumes that light partons in the incoming protons collide and produce pairs 
of heavy quarks QQ (Q denotes c or b quark) predominantly via the hard scattering 
processes of parton fusion qq — > QQ or gg — > QQ corresponding to the lowest-order 
(leading) graphs shown in Fig. El The next-to-leading order graphs of flavor excitation 
and gluon splitting mechanisms gain in importance only as the c. m. energy is considerably 
increased. 



Figure 2: Feynman diagrams for direct production of heavy quarks (Q is c or b quarks) via 
parton fusion mechanisms 

Therefore, we use PYTHIA 6.158 with the option MSEL = 5 to generate the bb events 
and MSEL = 4 for production of cc pairs (MSEL is a steering parameter in the PYSUBS 
common block). In this regime the simulation of heavy flavor production is performed 
only via the parton fusion mechanism with massive matrix elements for quark generation. 
Each event contains at least one QQ pair. The Gluck-Reya-Vogt (GRV94L) leading order 
proton parton distribution set and SLAC (Peterson) fragmentation function have been 
used, which are available when setting MSTP(51)=4 (by default) and MSTJ(11)=3, respec- 
tively. By default, the mechanism of neutral B- meson oscillations is switched on, which 
corresponds to setting MSTJ(26)=2 with oscillation parameters PARJ(76) = Xd = 0.7 
and PARJ(77) = x s = 20. 
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Let us distinguish three types of leptons originating from semileptonic weak transitions 
of heavy quarks - direct and indirect leptons from bb events, and direct leptons from cc 
events - according to the classification illustrated in Fig. 01 In Fig. El only specific decays 
of heavy mesons are shown but a full sample of events generated by PYTHIA contains 
also leptons originating from decays of heavy quark in baryons as well as decays of other 
particles. 
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Figure 3: Classification of leptons induced by semileptonic decays of heavy quarks: a) direct 
leptons from b or b due to transitions b — > l~X or b — > l' + X'; b) indirect (cascade) leptons from 
b or b due to transitions b — ► l~ Xc{-^> l' + X') or b — > l + Xc(l'~ X') c) direct leptons from c or c 
due to transitions c — > l + X or c — > l'~X' 
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Table V. Fractions of events with charged leptons without and with applying a cut on 
lepton transverse momentum px for various types of event selection: a) events with one 
or more charged leptons; b) ji ± e T pairs; c) like-sign lepton pairs / ± / ± 



Event 


Values 


No pt cut 


p T > 1 GeV 


selection 










a bb l°bb 




26 


a) 


u cc 1 u cc 


0.34 


0.012 




(1>1) i (2>1) 
cc l u bb 


540 


45 




<b Kf 


0.089 


8.9 x 10~ 3 


b) 


" cc 1 w cc 


0.019 


2.8 x 10~ 5 




cc l u bb 


210 


3.2 




bb 1 bb 


0.096 


6 x 10~ 3 


c) 


w cc 1 w cc 


0.0014 


6 x 10" 8 




(lift), (lift) 

cc l u bb 


15 


0.01 



The rates of various types of bb and cc events with charged leptons in final state, 
estimated using PYTHIA, are shown in Tabled The ratios of cross sections of cc and bb 
leptonic events were calculated assuming a ratio of total cross sections a^f / a^f » 1000. 
The leptons in the final state of bb and cc events can arise from both decays of heavy 
mesons or baryons and other sources including decays of kaons and pions. Tabled shows 
that the cc background dominates not only in the case of e/x-pair selection but even for 
events with like-sign lepton pairs if no cut on lepton transverse momentum px has 
been applied. A significant relative decrease of cc background in the case of pair 
selection as compared with e/i pairs is caused by a fact that the contribution of doubly 
semileptonic decays of bb to Z*/* sample arises as a combination of direct and indirect 
leptons, 

b->rx & b->l + Xc(->l'-X') or b^rXc{^l + X) & l->i'-x', (1) 

and combination of direct leptons due to B° — B° oscillations, 

b — > rX & b ^ b ^ l~X or b->b->l'-X' & b->l'-X', (2) 

while the contribution of cc events to the l^l^ sample arises mainly as a combinatorial 
background involving decays of kaons and pions. 

Due to simple kinematical arguments the lepton transverse momentum px is larger 
for leptonic decays of heavier particles. Therefore, one can expect that the px cut can 
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significantly reduce the combinatorial cc background to the signal of doubly semileptonic 
decays of 66, especially, in the case of selection of like-sign lepton pairs. Table [T] shows 
such a reduction of the cc background and its complete elimination in the case of l ± l ± - 
pairs selection after applying the requirement px > 1 GeV. Thus, the selection of like- 
sign lepton pairs with the additional pt cut can be efficiently used to trigger the signal of 
doubly semileptonic decays induced by decays of bb. Therefore, to study the asymmetry of 
lepton pairs induced by B — B production asymmetry, we restrict ourself to consideration 
only of the / ± / ± pairs. 

Let us define the overall charge asymmetry of like-sign lepton pairs as 



For Monte Carlo studies, a sample of more then 2 x 10 6 bb events with Z ± / ± pairs in final 
state with px > 1 GeV has been generated by PYTHIA. The prehistory of each particle 
from an event can be derived from the information stored in PYJETS common block, and 
the origin of each lepton can be traced back to the parton level, i.e. to the decay of heavy 
quark inducing the lepton. In this way, in accordance with the classification of leptons 
given in Fig. each event can be characterized by a code number "UK" , where "I" is 
the number of direct leptons from 66, "J" is the number of indirect leptons from 66, and 
"K" is the number of leptons from other quarks. 

There are two types of events dominating in the generated sample with 66-induced 
pairs corresponding with the following code numbers: 

• "110" - with one direct and one indirect lepton originated from doubly semileptonic 
decays (JTJ) of 66; 

• "200" - contribution of B° - B° oscillations 0. 

The full Monte Carlo sample also contains a fraction of events (about 20% in total) with 
the code numbers "210", "120", "020", "101" and "011". The event statistics and results 
on the estimates of overall asymmetry (jSJ) for the full sample of 66-induced like-sign lepton 
pairs and the dominating fraction of "110" and "200" events are shown in Tabled for 
various pt cuts. The fraction of the "110" events is decreased while the fraction of the 
B° — B oscillation contribution "200" is increased when increasing the pt cut. 

The main source of 66 events contain decays of I?-mesons while the contribution 
of 6-baryon decays does not exceed 1%. Therefore the contribution of the 6-baryon 
production asymmetry to the total charge asymmetry of lepton pairs is small, and the 
latter one reflects mainly the B — ^-production asymmetry. If we define the overall 
B — 5-production asymmetry as 



A 



N(l+l+) - N(l~r) 

N(i+i+) + nJFF) ' 



(3) 



A{B q ) 



N(B^) - N(B ) 
N(B<to)) + N(B m ) ' 



(4) 
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Table 2: Total Monte Carlo statistics, fraction of "110" and "200" events, and overall 
charge asymmetries (J3J) for 66-induced like-sign lepton pairs for various pt cuts 



Pt cut 


Number 


Fraction (%) 


Overall asymmetry (%) 


(GeV) 


of events 


"110" 


"200" 


Total 


"110" 


"200" 


1 


3 251 390 


43 


41 


2.00 ±0.06 


-0.99 ±0.09 


5.79 ±0.11 


2 


128 987 


27 


70 


3.55 ±0.27 


-1.92 ±0.63 


5.84 ±0.44 


3 


5 128 


19 


78 


2.22 ± 1.4 


2.8 ±3.6 


0.5 ±1.7 



the simulation of S-meson production in pp- interactions at the HERA-B energy gives the 
following estimates: 

A(B° d ) » -0.3% , A{B ± ) « 4.2% , A(B° S ) « -10.2% . 

The lepton charge asymmetries in the restricted phase space regions, i.e. as a functions 
of lepton transverse momentum px, Feynman variable xf and rapidity y of Z ± / =b pairs, 
are shown in Figs. El El and El Within errors due to limited Monte Carlo statistics, there 
is no noticeable px dependence in the lepton charge asymmetry of the l^l^ events with 
the code number "200" , caused by B° — B° oscillations. The considerable pr dependence 
in the lepton charge asymmetry of total sample of events reflects the strong dependencies 
of the contribution of "110" events on this kinematical parameter. The shapes of the 
corresponding histograms for "110" events in Figs. EJ El and El reproduce qualitatively 
the pt, if and y dependences of B — B production asymmetries. 




Figure 4: Dependence of charge asymmetry of (a) B-mesons and (b) 66-induced Z ± Z ± pairs on 
the lepton transverse momentum px 
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Figure 5: Dependence of charge asymmetry of (a) B-mesons and (b) 66-induced l^l^ events on 
the Feynman variable xf 




Figure 6: Dependence of charge asymmetry of (a) B-mesons and (b) 66-induced 11 events on 
the rapidity y 

To ensure that the observed effect of the charge asymmetry of the like-sign lepton pairs 
is really a manifestation of the B — B production asymmetry, we have repeated for the 
case of heavy-quark production in pp collisions a similar analysis of like-sign lepton pairs 
with the same selection criteria. Because of beam remnant symmetry in pp interactions, 
there is no B — B production asymmetry, and no charge asymmetry of lepton pairs was 
observed. 



Summarizing, we have found that at the proton energy E = 920 GeV, HERA-B 
provides a unique opportunity to study the B — B production asymmetry, caused at the 
fragmentation level by effects of asymmetric beam remnants for b and b quarks, by direct 
measurements of charge asymmetry in the production of like-sign lepton pairs without 
reconstruction of 5-meson decays. Our estimates shows that during a one year data 
taking run the statistics on Z*/ events at HERA-B would be large enough to provide a 
statistical error on charge asymmetry measurements at a level of a few percent even in 
the bins at the edges of the histograms in Fig. |S] and El corresponding to the borders of 
phase space in terms of variable xp and y where the asymmetry reach maximal values. 
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studies, useful discussions and careful reading of the manuscript. We are also grateful to 
T. Sjostrand and J. Carvalho for stimulating information and discussions. 
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